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ABSTRACT 



The t&val Postgraduate School computer simulation model QRAD was 
used to study sputtering from clean and nitrogen reacted tungsten, and 
from nitrogen reacted molybdenum, for (001) surfaces bombarded by 
normal incidence argon ions in the energy range from 0.5 keV to 
3.0 keV. The simulation results are compared to experimental data 
obtained by Winters, who concludes, that the difference in the cross 
sections for the sputtering of nitrogen from W(001) and Mo(001) is 
dominated by the adsorbate-substrate mass difference since the other 
important parameters are similar. The simulation however shows little 
dependence on the substrate mass, but a strong dependence of the 
sputtering cross sections on the distance of the adatoms above the 
substrate. The nitrogen position is required to be higher on tungsten 
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I. INTRODUCTION 



A. OVERVIEW 

Sputtering is the removal of surface atoms due to energetic 
particle bombardment. It is caused by collisions between incoming 
particles and atoms in the selvage, i.e., the near surface of solid 
material. First discovered by Grove [1] in 1853 and Faraday in 1854, 
it took about 50 years until the physical process involved vas 
recognized and another 50 years until a quantitative theory began to 
develop. Grove [1] noticed the disintegration of cathodes in glow- 
discharge tubes and observed, that the cathode material was deposited 
on the surrounding glass walls. He called this process "cathode 
sputtering". 50 years later, Goldstein [2] presented evidence that the 
sputtering effect was caused by positive ions of the discharge hitting 
the cathode. 

Stark [3,4] introduced a so-called 'hot-spot' model of sputtering 
(now known as 'evaporation theory') and later a collision theory, which 
viewed sputtering as a sequence of binary collisions initiated by the 
bombarding ion. Stark found that at low energies the yield, Y, 
increases with increasing ion energy, and interpreted this result on 
the assumption that an increasing amount of energy is transferred to 
the target atoms near the surface. At higher energies, increasing 
penetration depth of the ion causes a relatively smaller sputtering 
yield and thus a plateau in the yield curve. The reason for this is 
a decrease in the "size" of the ion as its energy increases. Stark 
considered his hot-spot model, the evaporation of target material from 
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a microscopically small region with high local temperature, and his 
collision theory as two different views of the same process. 

Von Hippel and Blechschmidt [5], who developed a quantitative form of 
Stark's hot spot model, made the first attempt to formulate a 
sputtering theory on the basis of local heating. Subsequent investi- 
gators found that this model is inadquate. Wehner [6] showed that 
local evaporation alone couldn't explain sputtering phenomena from 
single crystals. In 1934 Lamar and Compton [7] suggested that binary 
collision processes dominate in light-ion sputtering and local 
evaporation is dominant in heavy-ion sputtering. Their ideas led to 
the 'thermal spike' concept, in which a long lived high temperature 
volume presists in the target after the collision cascade is complete. 
Wehner' s observations revived the interest in the collision theory. 
Keywell [8] applied neutron transport theory to Stark's collision 
model. Due to his and subsequent calculations by Harrison [9] , 
probability concepts and collision cross sections were introduced as 
important parts of the sputtering theory. 



B. PHYSICAL MEANING 

Sputtering mechanisms operate in the near surface region of the 
target, and are caused by a collision cascade. The processes involved 
are simular to those causing radiation damage in the bulk material of a 
solid. Sputtering happens far from thermal equilibrium, which means 
that it isn't to be confused with evaporation. Generally an incident 
particle collides with atoms of the bulk material and transfers energy. 
If the transffered energy is larger than the cohesive energy at the 
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lattice site a primary recoil atom is created. This recoil atom 
collides with other target atoms and distributes the energy into a 
collision cascade. The condition for an atom to be ejected is that the 
kinetic energy associated with the normal component of its velocity is 
larger than the potential energy binding it to the surface. The 
sputtering yield is defined as the average number of target atoms 
ejected per incident particle. 

C. LATER WORK 

Although considerable discrepancies exist in the experimental 
results produced by various scientific groups, standards have developed 
which define the basic requirements for reproducible and reliable 
experiments. This development is shown in reviews of Behrisch [10] , 
Kaminsky [11], Carter and Colligen [12], and Pleshivtsev [13]. Over 
the last 30 years theoretical efforts have been increasingly successful 
in explaining the main features of sputtering. 

The analytic theoretical work in this area was done by Sigmund 
[14]. He presented a systematic theory of sputtering on the basis of a 
minimum of assumptions. Basic assumptions concern the characteristics 
of a single collision (differential cross section, inelastic stopping) , 
the structure of the medium (random or crystalline), the structure and 
binding forces of the surface, and some factors that depend on specific 
experimental arrangements like beam-target geometry, high or low ion 
energy, etc. The yield calculation consists of the following steps: 

1) to determine the amount of energy deposited by energetic 
particles (ion and recoil atoms) near the surface; 
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2) to convert this energy into a number of low-energy recoil atoms; 

3) to determine how many of these atoms come to the surface; 

4) to select those atoms that have enough energy to overcome the 
surface binding forces. 

Considering, for example, the backs puttering yield from a plane 
surface, the outward current of the target atoms through the surface 
plane at x = 0 is found to be: 

T dE d 2 fi 

J (E ,fl )dE d Z ft = TF-CE, 0, 0) — — |cos 0 | - 

0000 ° E IdE /dxl ° 4H 

o o 

where E is the energy of the ion impinging at an angle 0 onto a plane 
surface, F^ is the deposited energy density, ¥ is the number of 
projectiles per unit time, is the unit solid angle, and T m is a 
parameter which depends on Y and the scattering cross section. The 
sputtering yield is found by integrating this equation over E Q and 
0 O and dividing the sputtered current by the incident current 
This leads to 

Y = A F d (E, 0, 0) 

with 

A = **m / dE o / d (cos 0 ) | cos 0 | P (E , 0 ) 

T E IdE /dxl ° o o o 

o o 



where P(E 0 , 0 O ) is the probability that an atom will escape from 
the surface. The simplest model for the surface binding of a metal is 
based on a planar surface barrier U Q . In this case, 

1 for E q cos 2 0 o > U Q 



PCE 0 >©o) * 



0 for E rt cos^0 rt < U r 
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which leads to a simple form: 




4H 2 N C U 
o o 



with 



C = 
o 




2 

a BM» 



where A 0 = 24 and ag^ = 0.219 A. 

D. APPLICATIONS OF SPUTTERING 

For a long time sputtering was regarded as an undesired effect 
which destroyes the cathodes and grids in discharge tubes. In fusion 
research, plasma contamination by metal atoms is still one of the major 
problems. Sputtering also may cause the distruction of diaphragms and 
targets in accelerators and high-voltage electron microscopes. 

On the other hand, sputtering has become an indispensible process 
in modern technology. It allows the removal of very tightly bound 
surface layers on a nearly atomic scale. Sputtering is used in 
sputter-ion sources, for obtaining reproducibly clean surfaces and, for 
example, in micro-machining and depth profiling of thin films. 

Sputtering also can be used to obtain information about 
processes occuring on a atomic scale. Many scientists, including 
Kingdon and Langmuir (15), Winters (16), and Taglauer and Heiland (17), 
studied the sputtering of chemisorbed gas to investigate mechanisms 
which lead to the emission of multi-component materials. This thesis 
uses a computer simulation to varify one assumption made by Winters. 
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E. COMPUTER SIMULATIONS 



Due to their high operating speed, computers are attractive for 
calculations concerning sputtering events. A computer simulation can 
never replace an experiment, but by matching the results, a simulation 
can help to understand and identify sputtering mechanisms, and to 
interpret experimental results. A computer simulation is not a theory, 
it is a mathematical tool used to test the fundamentals of a theory 
[ 19 ]. 

Computer simulations can be divided in two categories: time-step 

and event-store models. Event-store programs move from event to event 
skipping the intervening time. For sputtering, they are based on 
binary collision theory. These fast, complex programs are useful when 
the physical model is well understood. The program used at the thval 
Postgradua te School is a time-step model based on simultaneous multiple 
interactions (Ml). MI simulations of physical systems usually require 
the calculation of about 100 individual ion trajectories (cascades) , 
while binary collision simulations use several thousands. This, and 
the fact that in MI simulations mechanisms are easy to trace, makes the 
time-step model attractive. 

This thesis is based on QRADS, a refined version of QRAD, the 
standard program in use at the Naval Postgraduate School for some time, 
which is a full-lattice simulation computer program, designed to model 
the dissipation of an incident ion's momentum in a single crystal 
target, using classical mechanics. The program has the ability to 
generate different surfaces of several crystal structures. To obtain 
reacted surfaces, a special subroutine places adatoms on the surface. 
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After the lattice is build, the potential functions are calculated. 

The impotant mathematical manipulations, concerning the trajectory, are 
done by the subroutine INTEG (integrator). This subroutine also tests 
for the termination of the trajectory. Finally, subroutine LEGIT 
determines which atoms were actually ejected. 



II. OBJECTIVES 



Some experimental results exist for sputtering from 
polycrystalline BCC structures [19] but only a few for clean single 
crystal surfaces or chemically reacted single crystal surfaces [19]. 
Winters [18] investigated the sputtering cross sections of tungsten and 
molybdenum bombarded by argon and xenon ions in the energy range from 
0.3 keV up to 5.0 keV. Figure 1 shows the cross sections as a function 
of incident ion energy. Winters' experiments show that the sputtering 
yield of nitrogen tends to increase as the atomic weight of the 
surrounding atoms increases. He furthermore concludes, that the 
difference in the cross sections for the sputtering on nitrogen from 
W(001) and Mo(001) is dominated by the absorbate-substra te mass 
difference since the other important parameters are similar. 

The primary objective of this thesis is to study sputtering from 
nitrogen reacted (001) tungsten and molybdenum surfaces using a 
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computer simulation model and to compare the results with the exper- 
imental data obtained by Winters. 

The incident ion will be Ar + in the energy-range from 0. 5 keV to 
3.0 keV. 

The simulation includes: 

1) Bombardment of a clean (001) tungsten surface (used for 
comparison) . 

2) Bombardment of a nitrogen reacted (001) tungsten surface. 

3) Bombardment of a nitrogen reacted (001) tungsten surface where 
the atomic mass of tungsten is replaced by the atomic mass of 
molybdenum. 

4) Bombardment of a nitrogen reacted (001) molybdenum surface, where 
the nitrogen is closer to the target surface plane than in the 
cases listed above. (The nitrogen is assumed to be in a 4-fold 
symmetry position for all reacted surfaces.) 



III. THE SIMULATION 



A. THE TARGET 

1. Properties of Materials 

a) Tungsten (W, Wolfram): atomic mass 183.85, lattice parameter 

3.165 A, crystal structure BCC, nearest-neighbor distance 
2.74 A, cohesive energy 8.9 eV, 
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b) Nitrogen (N): atomic mass 14.0, binding energy to the tungsten 

substrate 6.5 eV, 

c) Molybdenum (Mo): atomic mass 94.94, lattice parameter 3.147 A, 

crystal structure BCC, nearest-neighbour distance 2.73 A, 
cohensive energy 6.82 eV. 

2. Location of the Adatoms 

The location of nitrogen atoms on a tungsten or molybdenum 
surface is not known. Clavenna [21] and Adams [22] suggest that the 
nitrogen is single-bounded in a 4-fold position above an atom in the 
second layer of the tungsten substrate. (The 4-fold position is also 
assumed by Winters.) Figures 2a, b show side and top view of the 
reacted surface where the radius of the nitrogen is assumed to be 0.59 
A [22]. (The smaller circles are the nitrogen atoms.) Griffiths 
et. al. [23], agree with this position at a fractional coverage (Of) 
of 0.5 but mention, that at this coverage the nitrogen is unstable. 

They recommend study at the most stable coverage of Of = 0.4. This 
requires formation of islands of nitrogen atoms on the surface, (see 
Figs. 3a, b). Harrison et. al. [24], propose that the nitrogen might be 
below the first layer of the substrate in an interstitial site. In 
this case some second layer atoms would be pushed down into the 
substrate to set equal distances to all nea rest-neighbors. For this 
simulation the nitrogen is assumed to be in the 4-fold position for 
both substrate materials. For the nitrogen reacted tungsten surface 
the location above the substrate is chosen so that the distance to all 
five nea rest-neighbors is equal. The distance above the surface is 
0.504 LU (where 1. 0 LU = 1. 58 A = 0. 5 lattice parameters, see Fig. 4). 



14 



(Lattice units were introduced in the program so that rescaling to a 
new target material can be done without major reprogramming and so that 
results from different substrates can be compared directly.) Runs also 
were made with the tungsten mass replaced by the molybdenum mass to 
examine the mass effect proposed by Winters [18] , discussed in chapter 
2. On the nitrogen reacted molybdenum surface the nitrogen is placed 
closer to the surface layer (see Fig. 4). The distance above the 
surface is 0.243 LU (where 1.0 LU = 1.57 A = 0.5 lattice parametrs). 

3. The Target Size 

The inability of the lattice of target atoms to contain the 
entire cascade generated by the incident ion is called ' a failure in 
containment'. The size of the target determines the containment of the 
trajectory. If an increase in target size doesn't change the results 
of the computations, the trajectory is effectively contained for the 
purposes of the computations. (Unnecessarily large targets waste 
computer resources. ) 

The minimum target size for the ion energy of 0.5 keV vas found 
to be 19x8x19, where the numbers represent the numbers of planes in the 
x,y and z directions. 40 nitrogen atoms (equivalent to a fractional 
converage of 0f = 0.4) are placed as described above (see Fig. 5; 
solid lines mark the smaller target). The minimum target size for ion 
energies between 1.0 keV and 3.0 keV v&s found to be 23x8x23 with a 
fractional coverage of Qf = 0.416, equivalent to 60 adatoms. The 
difference in coverage does not significantly affect the result. 

Winters [18] , determined the coverage in his experiments by Auger 
spectroscopy, and found that the Auger intensity is a nonlinear 
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Thus all 



function of coverage above 9 = 6 x 10 3 ^ atoms per cm 2 , 
his measurements were taken below this value. 

The coverage chosen for the simulation is equivalent to 4.9 x lO 3 ^ 
atoms per cm 2 . 

B. POTENTIAL FUNCTIONS 
1. In General 

A primary uncertainty in sputtering models and theory is the 
form of the interatomic potential functions. No single potential 
function has yet been deduced which fits every sputtering experiment. 
The functions used for this simulation are Morse, and Moliere [25] 
joined with a cubic spline [24]. The Morse potential function is of 
the form: 

V(r) = D e exp[-2a(r-r Q ) ] - 2D e exp[-a(r-r Q ) ] 
where D e is the well depth, a is a scale factor of the potential 
function, and r Q is the equilibrium separation. For r > r Q the 
Morse potential is attractive and for r < r Q repulsive. 

The Moliere potential function is of the form: 

V(r) = (Z x Z 2 e 2 /r) (0. 35 exp[-0.3r/a] + 0.55 exp[-l. 2r/a] 

+0.10 exp [-6. Or /a]) 

where 'a' is the Firsov screening length. The Moliere function is a 
pure repulsive potential function. 

The cubic spline is of the form: 

V(r) = C Q + C x r + C 2 r 2 + C 3 r 3 

where the coefficients are determined by the program. The cubic spline 
is necessary to combine Morse and Moliere functions into a single 
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composite potential. For bard collisions, at close approach in the 
keV-range, the Morse potential bends towards the Y-axis and thus cannot 
be physically correct, because all potential functions must approach 00 
as r + 0 like 1/r. The cubic spline is inserted where the Morse and 
the Moliere functions overlap, thus the combined function is a Moliere 
to the left and a Morse to the right of the spline. In general the 
spline shouldn't be longer than 0.05 LU to get smooth potential and 
force curves. 

2. Crystal Potential Functions 

The composite potential functions are used for the bulk 
materials as well as for the reacted surfaces. The Moliere function is 
a standard Moliere in each case. The Morse parameters have physical 
meaning: a, used as an adjustable parameter, is chosen so that the 

slope of the Morse matches that of the Moliere function. For the 
substrate materials D e is chosen to generate individual atom energies 
in the center of the target equal to the cohesive energy, i.e., 8.9 eV 
for tungsten and 6.82 eV for molybdenum. For the reacted surfaces D e 
has to match the binding energies between the nitrogen and the substate 
i.e., 6.5 eV in each case. Because the potential function is truncated 
after one nearest-neighbor distance, r Q is greater than the 
separation between atoms in the substrate materials. For the reacted 
surfaces, r Q is the nearest-neighbor distance. 

The cubic spline is bounded by the values and Rg. The 
table below shows the parameters for the crystals amd adatoms. 
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Table 1: Parameters for Crystal Potential Function. 





oCA* 1 ) 


r e (A) 


D e (eV) 


r a ( a ) 


R b (A) 


see Figs. 


w -w 


1.20 


2.849 


1.335 


1.83 


1.90 


6,10 


W -N 


1.65 


2.380 


1.300 


1.14 


1.14 


7,10 


Mo -Mo 


1.50 


2.800 


0.997 


1.24 


1.30 


8,10 


Mo-N 


2.20 


1.950 


1.612 


0.91 


0.97 


9.10 


Figures 


6, 7, 8, 


9 show 


the Morse and 


Moliere 


functions and Figure 



10 shows the combined functions generated by the program. 
3. Gas Phase Substrate Atom Potential Functions 



The gas phase potential functions are pure Morse functions 
because no hard collisions are expected. The parameters are listed 
below: 



Table 2: Parameters for Gas Phase Substrate Atom Potential Functions. 





<*(A _1 ) 


r e (A) E 


> e (eV) 


see Figs. 


W -W 


1.20 


2.74 


8.900 


11 


W -N 


1.65 


1.96 


6.721 


11 


Mo -Mo 


1.50 


2.73 


6.820 


12 


Mo-N 


2.20 


1.96 


5.100 


12 


4. 


Nitrogen- 


•Nitrogen 


Potential 


Functions 



The N-N potential function in the gas phase as well as in the 
crystal phase is a pure Morse functions with: 



= 2.7 


A-l 


e = 1.098 


A 


i e = 7.373 


eV (crystal phase) 


i e = 9.834 


eV (gas phase) 



see Figures 11, 12 for the gas phase potential functions. 

5. Argon-Tungsten/ Molybdenum and Argon- Nitrogen Potential 



Functions 

Because it's not likely that argon ions react with nitrogen 
tungsten or molybdenum, these potential functions have the standard 
Moliere form, (see Fig. 13). 
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C. THE IMPACT AREA 



The impact area for each simulation is a rectangle containing a 
grid of 104 impact points, (see Fig. 14), where the length of the sides 
is 3.165A = 2.0 LU for tungsten and 3.147 A = 2.0 LU for molybdenum. 

The origin of the rectangle is placed at the following locations: 



Table 3: Origns of the Impact Area 



target size 



locations 



23x8x23 

19x8x19 

23x8x23 

19x8x19 



N-reacted 

N-reacted 

clean 

clean 



(10.10) (10,12) Fig. 5 

(08.08) (08,06) Fig. 5 

(11.11) Fig. 5 

(09.09) Fig. 5 



On the Nitrogen-reacted surface, two impact point sets are neces- 
sary to assure complete converage. 



IV. DISCUSSION AND RESULTS 



All output vas analyzed by a separate analyzer program, AN83, which 
provides a numerical analysis, and by the program PLOTAN, providing a 
graphical analysis. The latter is discussed first. 



A. GRAPHICAL ANALYSIS 
1. Spot Pattern 

The spot pattern plots show the polar and azimuthal angular 
distribution of the ejected atoms, (see Figs. 15, 16, 17). The 
distances of the crosses form the center is a measure of the polar 
angle of the ejected atoms, i.e., a cross far out corresponds to a 
large ejection angle measured from the surface normal. The spot 
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pattern obtained for the substrate materials (W, Ho, Mo-mass) show a 
four-fold symmetry, where ’Mo-mass' means, that in the program the 
tungsten mass v&s replaced by the molybdenum mass and all other 
parameters were unchanged. 

The spots at 45°, 135°, 255° and 315° point in the direction of 
the nearest-neighbors. Figure 15 shows the spot pattern for molybdenum 
atoms ejected from the substrate after bombarding a nitrogen-reacted 
molybdenum surface by 2.0 keV argon ions. This figure is typical of 
all spot pattern obtained from substrate materials. 

The spot pattern obtained from the adatoms show no symmetry. 
Figures 16, 17 are representative for ejected nitrogen atoms. Figure 
16, the case where the nitrogen is placed further above the surface, 
shows a heavier concentration of ejected atoms at low angles than 
Figure 17. 

2. Energy Distribution 

Figure 18 shows typical energy distributions for atoms 
sputtered from the substrate in the ion energy range from 0.5 keV to 
3.0 keV. The envelopes peak near 7.0 eV. 

The energy distribution of the sputtered nitrogen atoms depends 
on the location of the adatom. If the nitrogen is close to the surface 
layer, the distribution peaks at about 2.5 eV/atom, (see Fig. 19). If 
the nitrogen is further away, as in the case of the tungsten or 
molybdenum-mass substrate, the distribution is much broader, and the 
maximum is roughly 15 eV/atom, (see Fig. 20). 
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3. Energy Distribution >100 Fsec 



The energy distribution for particles ejected from the 
substrate after 100 fsec (1.0 fsec = 10“^ sec) is similar for all 
materials, but with the difference that the fraction for molybdenum or 
molybdenum-mass is less than for tungsten, (see Figs. 21, 22, 23). 

More low energy atoms are ejected after 100 fsec in the case of the 
heavier substrate material. This could indicate a mass effect. 

Energy distributions for the ejected nitrogen atoms are shown 
in Figures 24 and 25. If the nitrogen is in the lower position there 
were no atoms ejected with energy greater than 7.0 eV. This behavior, 
is seen throughout the whole energy range, (0.5 keV to 3.0 keV). 
Although the total number of atoms is low, (see Fig. 25), and therefore 
the statistics are poor, the difference between Fig. 24 and Fig. 25 may 
be indicating a distance effect. 

4. Energy Versus Time 

Figures 26, 27, 28 and 29 show that no high energy particles 
are ejected later than 100 fsec. The spectrum for substrate materials 
(Fig. 26) is wider than that for the ejected nitrogen atoms. 

Figure 28 shows that hardly any nitrogen atoms are ejected 
after 100 fsec. The difference between Fig. 27 and 29 and Fig. 28 on 
the other hand shows again, that there may be a distance effect. 

5. Ejection Time Distribution 

The ejection time distribution is essentially the same for all 
substrate materials. For high energies the distribution is lower and 
broader, (see Figs. 30 and31). The distribution peaks at about 80 
fsec. 



21 



6. Atoms per Single Ion 



The atoms per single ion diagrams are similar throughout all 
runs. The spectrum broadens with increasing energy, (see Figs. 38, 39, 
40 and 41) . 

7. Yield per Single Impact Point 

The yield per single impact point diagrams are also similar for 
all runs, (see Figs. 42 and 43), where the open circles represent no 
yield. The diagrams are combined for the locations (10, 10) and (10, 
12), i.e., the upper left corner of each rectangular box represents the 
location (10, 10). Comparison of the centers of the combined squares 
shows, that whenever the ions hit a substrate atom the yields are 
larger. The center of each right square is occupied by a nitrogen 
atom. The nitrogen acts like a 'buffer'. The open circles are caused 
by channeling. Figures 42 and 43 show that with decreasing ion energy 
the areas where channeling occurs decrease. The reason for this is 
that for decreasing energies the cross sections increase because the 
impact energy appears in the denominator. 

8. Yield per Surface Atom 

Figures 44 , 45 show the yield obtained from the surface layers 
for both target sizes. Almost no atoms are ejected from the edges of 
the crystal. This indicates that the crystals are large enough to 
contain the trajectories. 

B. NUMERICAL ANALYSIS 

The numerical analysis done by the analyzer program AN 83 provides 
the following nitrogen yields, (see also Fig. 46): 
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Table 4: Nitrogen Yields from Simulations. 



substrate : 


W(001) 


Mo-mass(OOl) 


Mo (001) 


energy 
in keV 


yields 


yields 


yields 


0.5 


2.36 


2.04 


0.60 


1.0 


2.25 


2.29 


0.68 


2.0 


1.80 


2.19 


0.51 


3.0 


1.37 


1.87 


0.51 



Experimental yields in the ion energy range from 0.5 keV to 3.0 keV 
are known only for polycrystalline tungsten and molybdenum surfaces 
[26]. Figure 47 shows the experimental metal yields obtained from 
polycrystalline tungsten and molybdenum as well as the yield from clean 
W(001) obtained from the simulation. The yields are listed below: 

Table 5: Yields from Clean Metal Surfaces. 

Experiment Simulation 



energy 
in keV 


W(poly) 


Mo( poly) 


W(001) 


0.5 


0.58 


0.57 


2.31 


1.0 


0.88 


0.80 


2.94 


2.0 


1.30 


1.10 


2.40 


3.0 


1.50 


1.30 


2.18 



Some data are known for 5.0 keV Ar + bombardment of molybdenum 
[20] at normal incidence onto the (001) planes and onto polycrystalline 
molybdenum. The yield from polycrystalline Mo is found to be larger 
than from Mo(001). Thus Roosendaal [20] concludes that for a BCC 
crystal the same result is found as for FCC structures: the sputter- 

ing yield for less closed packed planes is less than the yield from 
poly crystalline targets, i.e., Y(Poly) > Y(001). The results of this 
simulation show a yield of about twice the yield from polycrystalline 
tungsten. Winters [16] and Wehner [27] who bombarded W(001) surfaces 
in the low energy range with various ion types obtained the results 
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shown below 



Table 6: Experimental Yields from Clean Metal Surface for Ion 

Energy of 0.5 keV. 





Neon 


Argon 


Krypton 


Xenon 


Winters : 


0.25 


— 


0.91 


1.21 


Wehner : 


0.27 


— 


0.92 


1.02 



The yield is increasing with increasing ion mass and should be 
about 0.6 for 0.5 keV Ar + ions. The result of this simulation is 
higher than this value by a factor of about 4. This shows that the 
(Ar-W) potential (a standard Moliere) is too hard in this energy 
range. 

Winters [18] shows his results as cross section versus ion energy. 



(see Fig. 


1) . The nitrogen cross 


sections are: 


Ta ble 7: 


Experimental Nitrogen Cross Sections. 

cross sections x 10”^^, 


energy 
in kev 


substrate W(001) 


Mo(001) 


0.5 


0.91 


0.65 


1.0 


0.98 


0.72 


2.0 


1.00 


0.75 


3.0 


0.94 


0.73 



Winters assumes that the nitrogen sputtering rate is given by 

~ Rjj = dO^j/ dt = 

where ©jj is the nitrogen atom concentration (atoms/cm^), Opj 
is the sputtering cross section (cm^), v + the ion flux (ions/cm*" 
s) , and t the time in seconds. The product tv + is the ion dose 
(fluence) while the product gives the sputtering yield in 

nitrogen atoms/ ion. 
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Using this relation (Y = the cross section obtained from 



the simulation are: 



Table 8: Simulation Nitrogen Cross Sections. 





cross sections x 


10" 15 , (cm 2 ) 




energy 
in keV 


substrate: W(001) 


Mo-mass(OOl) 


Mo (001) 


0.5 


4.77 


4.12 


1.20 


1.0 


4. 55 


4.63 


1.36 


2.0 


3.64 


4.42 


1.02 


3.0 


2.77 


3.78 


1.02 



Figures 48 and 49 show Winters cross sections in comparison to the 
cross sections obtained from the simulation. In the case where the 
nitrogen is in the lower position (Mo) (Fig. 48) , the simulation cross 
sections are in a better agreement with the experimental results. But 
even in this case the Simula ted cross sections are still too large. 

This suggests that the (Ar-N) potential (also a standard Moliere) is 
too hard. Winters' experiments showed that the sputtering yield of 
nitrogen tends to increases as the atomic weight of the substrate 
increases. He suggests that the difference in the cross section for 
the sputtering of nitrogen from W(001) and Mo(001) is dominated by the 
adsorbate-substrate mass difference since the other important 
parameters are similar. The simulation shows that the distance between 
the nitrogen and the substrate is more important than the mass 
difference. 

Winters furthermore states that the shape of the yield curves for 
the chemisorbed gas and the substrate material is different because the 
ratio a (gas)/o (substrate) decreases with increasing ion energy. 

Winters used a (substrate) cross sections obtained from polycrystalline 
W and Mo assuming a surface coverage of 9 = 10 + ^ atoms/cm“. The 
following ratios were calculated from the simulation: aN(W)/aW(clean) , 
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a (Mo-mass) /a W(clean) , oN(Mo)/oW(clean) . Using aW(clean) in the 
latter case seems reasonable, because comparision of yields from 
polycrystal- line W and Mo as well as the yields obtained from the 
substrates are not too different, i.e. , the yields from Mo(clean) 
shouldn't diviate much from W(clean). 0 for W(001) is 9.98 x 10 + ^ 
atoms/cm^. The ratios are shown below (see also Fig. 49). 

Table 9: Cross Section Ratios. 







Experiment 






Simulation 




oN(W) ' 


oN(Mo) 


oN(W) 




a N( Mo-mass) 


aN(Mo) 


energy o 


W(poly) 


oMo(poly) 


aW(clean) 


aW(clean) 


aW(clean) 


in keV 














0.5 


1.57 


1.14 


2.06 




1.78 


0.52 


1.0 


1.11 


0. 90 


1.54 




1.57 


0.46 


2.0 


0.77 


0.68 


1.52 




1.84 


0.43 


3.0 


0.63 


0.56 


1.27 




1.73 


0.47 


In comparison to 


Winters' results, 


Figure 


50 


shows the ratios for 


nitrogen 


in the 


lower position and 


Figure 


51 


for nitrogen 


in the higher 


position. 


(Mo-mass not shown.) The values 


obtained from 


the Simula- 



tion show the general trend observed by Winters for the W and Mo 
substrate, but indicate that the location of the high nitrogen is 
perferable for the tungsten substrate and low for the molybdenum 
substrate. Even with these choices, the difference between the 
nitrogen locations is too large. This is also seen by comparing 
aN(W)/aN(Mo) between Winters data and the data obtained from the 
simulation. The ratios are shown below, (see also Fig. 52): 



Table 10: Cross Section Ratios. 



energy 



Experiment 



Simulation 



in keV 


oN(W)/oN(Mo) 


0.5 


1.40 


1.0 


1.36 


2.0 


1.33 


3.0 


1.29 



aN(W)/aN(Mo-mass) aN(W)/aN(Mo) 



1.16 


3.98 


0.98 


3.35 


0.82 


3.57 


0.73 


2.72 
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As mentioned in section A7, with decreasing ion energy the cross 
sections increase and the areas where channeling occur decrease* This 
is also seen comparing the reflected argon yields, listed below: 



Table 11: Simulation Argon Yields. 





clean 


nitrogen 


reacted 




energy 
in keV 


substrate W(001) 


Mo-mass(OOl) 


W(001) 


Mo(001) 


0.5 


0.40 


0.30 


0.43 


0.24 


1.0 


0.27 


0.11 


0.22 


0.13 


2.0 


0.12 


0.07 


0.09 


0.05 


3.0 


0.07 


0.07 


0.05 


0.03 


Due to 


the increase in cross 


sections , more ions 


are reflected at 


low energies. Harrison (19) 


says that about 85% 


of the 


ejected atoms 



(bulk material) from clean targets come from the first layer. This 
simulation shows that at least 80% of the ejected atoms come from the 
first layer. The yield from deeper layers varies between 19% and 5%, 
where the yield from below the second layer is almost neglible. The 
yields from below the first one are shown below: 



Table 12: Yields from Second or Deeper Layers (Simulation), 

clean nitrogen reacted 



energy 
in keV 


substrate 


W(001) 


Mo-mass(OOl) 


W(001) 


Mo (001) 


0.5 




11% 


12% 


16% 


5% 


1.0 




10% 


14% 


17% 


8% 


2.0 




16% 


16% 


19% 


12% 


3.0 




17% 


18% 


17% 


16% 


Some of 


the ejected 


atoms, which 


have similar 


energy , 


velocity and 



position, combine above the surface, forming mul timers. Most of the 
multimers observed are dimers: (W-W, (W-N) , Mo-Mo) , (Mo-N). No (N-N) 
dimers were found. The numbers of dimers observed are listed below: 
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Table 13: Dimers from Clean Tungsten (Simultion) . 



energy 


W - 


in keV 




0.5 


40 


1.0 


53 


2.0 


18 


3.0 
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Table 14: 


Dimers from 


Nitrogen Reacted 


Surfaces 


(Simulation) 


energy 


W - N 


52 

1 

O 

S 


W-W 


Mo - Mo 


in keV 










0.5 


50 


9 


18 


18 


1.0 


65 


12 


23 


31 


2.0 


27 


15 


24 


15 


3.0 


14 


8 


19 


10 



The table below shows the number of linked atoms (taken from all 
mul timers) . 



Table 15: Linked Atoms (Simulations). 

energy substrate W Mo- mass Mo 



in keV 


(WW) 


(WN) 


(MoMo) 


(MoN) 


(MoMo) 


(MoN) 


0.5 


33 


71 


56 


85 


33 


18 


1.0 


41 


82 


74 


92 


54 


29 


2.0 


72 


66 


37 


64 


27 


18 


3.0 


42 


37 


42 


65 


29 


19 



In the case of the molybdenum substrate the number of Mo-N links 
exceeds the number of Mo-Mo links throughout the whole energy range. 
In the case of the tungsten substrate the number of linked W-W is 
higher for 3 keV and 2 keV. 
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V. SUMMARY AND RECOMMENDATIONS 



The results of the simulation show that the ion-substrate 
potential functions are chosen too hard. Further work should be done 
to reduce the yields for the substrate materials. The graphical 
analysis shows that there may be mass as well as distance effects, but 
the results of the simulation clearly show that the adsorbate-substrate 
distance is more important than the adsorbate-substrate mass 
difference, when the nitrogen sits in the higher position on the 
substrate with the higher atomic mass. After adjusting the substrate 
yields, futher work should be done to find the right distance between 
the adatoms and the substrate. 
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FIGURES 




F)C. I Sputtering eras sections for nitrogen from WMOOl and McV ITJOJ 
versus ton energy. The sputtering yield is the product of the eras section 
and the atomic density of adsorbed nitrogen. 



Fig. 1. Sputtering Cross Sections. 
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a) 



N W 




Side view 



b) 




Topview 

Fig. 2. Nitrogen reacted tungsten surface. 

Small circles represent the N-atoms. 
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' w* Vjl 'wL v 



for the C . 4 monolayer IC on ’V(CCI) 
surface: Large hatched circles, top-layer 



"’-atoms, small filled circles, F-atcr.s. 

a) Plan view, illustrating domain and 
boundary structure. 

b) Cross section through the domain 
along AAl 



Fig. 3. Island formation of a nitrogen reacted 
tungsten surface. 
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a) 



Sideview 




b) Topview 

A s 




Small circles represent nitrogen atoms, 

* 

Fig, 4* Location of nitrogen atoms on tungsten 
and molybdenum substrates. 
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Small circles represent nitrogen atoms. 
Solid lines mark smaller target. 

Numbers represent layers. 

Fig. Topview.. of ..reacted surface. 
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Fig. 7. Interatomic potential functions. 
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potential energy, ( eV ) 




O 

Separation, (A) 

Fig. 8. Interatomic potential functions. 
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potential energy, ( eV ) 



A 




Fig. 9 * Interatomic potential functions. 
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Fig. 10. Interatomic potential functions. 
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potential energy, ( eV ) 




0 1.0 2.0 , 3.0 4. 

Separation, (A) 

Fig. 13. Interatomic potential functions. 
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Fig. 1^. Impact area. 
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2.0 KEV M0(001) 'MO' 



A. 




SPOT PATTERN 

(obtained for Mo-substrate) 




46 





3.0 KEY W(00t) *N* 



SPOT PATTERN 




Fig. 16. Spot pattern for nitrogen atoms from a 
W(C01) surface. 
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3.0 KEV M0(001) 'N' 



SPOT PATTERN 




♦ 

* 

Fig. 17. Spot pattern for nitrogen atoms from a 
Mo (001) surface. 
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ENERGY (EV) ENERGY (EV) 

Fig. 18. Energy distribution for ejected Mo atoms from a reacted surface 
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0.0 6.0 10 0 16.0 20.0 26.0 

ENERGY (EV) 

Fig. 20. Energy distribution for ejected N atoms from reacted surfaces. 
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orgy distribution for atoms (substrate), ejected after 100 fsec 
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Fig. 22. Energy distribution for atoms, ejected after 100 
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Fig. 23. Energy distribution for atoms (substrate), ejected after 100 fsec 
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0.6 KEV MO(OOI) N* 



ENERGY DISTRIBUTION > 100. FSEC 




Fig. 25. Energy distribution for N acorns, ejected after 
TOO fsec. 
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Fig. 26. Energy versus time for W atoms ejected from a reacted surface 
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Fig* 27. Energy versus time for N atoms ejected from a reacted surface 
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tion time distribution for substrate atoms 
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31 • Ejection time distribution for substrate atoms 
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Fig. 32. Ejection time distribution for N 
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Fig. 33 * Ejection time distribution for N atoms 
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Fig. 34« Ejection time distribution for II atoms 
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Fig. 35. Ejection time distribution for N atoms 
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Fig. 37. Ejection time distribution for N atom 
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Fig. 33. ASI for ejected N atoms. 
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Fig. 39. ASI for ejected K atoms. 
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Fig. 'vO. ASI for ejected N atoms. 
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Fig, 41* ASI for ejected N atoms. 
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3.0 KEV M0{001) 'MO' 




3.0 KEV W(001) 'N' 




Fig. 42. Yield per single impact point. 
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1.0 KEV M0(001) 'MO' 




0.5 KEV MO(OOI) 'MO' 




Fig. 43« Yield per single impact point. 
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Fig. 45* Yield per surface atom. 
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Ion energy, (keV)' 

Fig. 47. Simulation yields from V/(001) and experimental yields 
from poly crystalline 'll and Ho. 
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Fig. ! i9* Sputtering cro i Lj 1 ;c Lions for adsorbed nitrogen 
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Vi(S. 50. Cx’oss sections ( n i Iro^xti ) / cross sections (substrata) 
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0.0 1.0 2.0 3.0 4.0 

Ion energy, (keV) 

Fig. 52. Cross section ratios for N from 17-substrate / N from 
lio-substrate. 
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